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A new method for the evaluation of the perturbative expansion of the QCD pressure 
is presented which is valid for all temperatures and quark chemical potentials in the 
^ • deconfined phase, and worked out up to and including order g 4 . This new approach 
unifies several distinct perturbative approaches to the equation of state, and agrees with 
. dimensional reduction, HDL and HTL resummation schemes, and the zero-temperature 
result in their respective ranges of validity. 

\Q : 

1. Introduction 

43 : 

Q_i! In a recent paper [ [T] Kajantie, Rebhan, Vuorinen, and myself have proposed a new 
method for calculating the QCD equation of state in a perturbative setup for all temper- 
atures T and quark chemical potentials /i. Previous perturbative approaches have been 
limited to certain regions only: when the temperature is large compared to the Debye 
screening mass T > m-Q (with m 2 ^ = g 2 (T 2 + fi 2 / n 2 ) for N c = 3 colors and Nf = 2 flavors), 
hard modes can be integrated out to leave an effective theory for the zero mode, and di- 
mensional reduction applies [I2j-[E]- I n the opposite case T < m D , a four-dimensional 
resummation in terms of hard dense loops (HDL) works, giving rise to anomalous con- 
tributions in entropy and specific heat of a non- Fermi liquid [E]- The latter approach 
smoothly connects to the well-known zero-temperature limit by Freedman and McLerran 
[ E] • Our new method covers all of these cases, and provides an independent check on 
those calculations. 

These parametrically distinct regions are suitably visualized in a QCD phase diagram 
in polar coordinates (Figs. ^El) • Due to asymptotic freedom, weak-coupling techniques 
apply at high temperatures or chemical potentials r »s \JT 2 + fi 2 . We label the abscissa 
in Figs. 121 and El by powers of the coupling g 2 , g 3 , ... and the ordinate by parametrically 
different regimes T ~ /i, g/i, g 2 fi, . . . which for larger r are displayed as horizontal lines 
(somewhat deviating from strict polar coordinates). In the following, we will describe the 
various existing perturbative approaches, and connect them by our new method. 
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2. Previous methods 



If the temperature T is larger than all other dynamical scales, one can apply effective 
field theory methods [ |2] and integrate out the degrees of freedom corresponding to non- 
zero Matsubara modes to obtain a simpler three-dimensional (dimensionally reduced) 
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Figure 1. QCD phase diagram in polar coordinates. 



effective theory, known as electrostatic QCD (EQCD). The effective theory approach can 
be applied a second time to remove also the massive longitudinal gluon from the theory, 
thus producing an effective three-dimensional pure Yang-Mills theory, which is called 
magnetostatic QCD (MQCD). Practically, the effective Lagrangians and effective masses 
m-E and couplings g^, Xe, gM, ■ ■ ■ can be obtained by matching a set of physical quantities 
between a generic Lagrangian and the full theory. The perturbative pressure can then be 
written as 

V = Pe(T, g) + PAi(m 2 E , g E , X E , •••) + Pg(9m, ■ ■), (1) 

where the first two terms correspond to the coefficients of the unit operator of EQCD 
and MQCD respectively, and the last term corresponds to the pressure of MQCD, which 
constitutes a fundamentally non-perturbative contribution to full QCD. Figure 121 shows 
at which orders the latter two terms start to contribute, by the lines labeled EQCD and 
MQCD. The plot shows for example that the non-perturbative contribution of MQCD 
which sets in at order g 6 for T > /i, only affects the pressure at order g 8 if the temperature 
is of the parametric order T ~ g/j,. So far the perturbative pressure has been obtained up 
to and including order g 6 lng and extended to finite chemical potential \x [H]. 

The other result depicted in Fig. |21 is the zero-temperature limit by Freedman and 
McLerran (FMcL) [ [BJ which is known up to and including order g A , with an error of 
order g 6 lng. At order g 4 dimensional reduction diverges as T — > and does not connect 
to FMcL, revealing obviously the breakdown of dimensional reduction. 

In the region T < g/i perturbation theory requires a different reorganization which can 
be most efficiently performed via HDL resummation (Fig. EJ). A small modification to 
HTL resummation increases the range of validity of this approach to larger temperatures. 
The interaction pressure Sp = p — p SB — (p — p S b)|t=o (Psb is the interaction-free Stefan- 
Boltzmann contribution) as well as entropy and specific heat show anomalous logarithmic 
enhancement of the leading contribution (a signature for non-Fermi-liquid behavior) and 
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Figure 2. Range of validity of dimensional 
reduction (DR) and the zero-temperature 
result (FMcL). 




Figure 3. Range of validity of HDL and 
HTL resummations, and the piy presented 
here. 



fractional powers of the temperature in the regime T < gfi (for Nf = 2, N c = 3) 

6p = ^ |ln^ - 0.7906 - 2.9401 ( ^) "* + 2.2312 + . . A. (2) 

The fractional powers are indicated by lines starting from T ~ go. and g 4 in Fig. |21 and 
form the correct continuation of the plasmon term oc Tm^ (EQCD line in Fig. |2J). 

3. The new approach 

Our new approach [|T] is motivated by the observation that only a subset of all diagrams 
that contribute to the pressure exhibit infrared divergent behavior in the limit T — > 0. 
The problematic set consists of the two-gluon reducible diagrams (2GR; vacuum bubble 
diagrams that can be separated by cutting two gluon lines). Up to order g 4 , the following 
classes of 2GR diagrams require resummation: 




where the black dots represent the full one-loop gluon polarization tensor. This is only 
known numerically, but the contribution of diagram (b) to the pressure up to 0(g 4 ) can 
be expressed analytically, and diagrams (c) and (d) would only contribute beyond that 
order. Only the contribution from diagram (a) has to be obtained by means of numerical 
integration (using the methods developed at large Nf [E]). The pressure through order 
g 4 (labeled by the roman numeral IV) valid for all T and is finally given by the sum 
of analytically calculable pieces and the ring diagrams (a) which have to be integrated 
numerically: p w = p an i + ^ a n fe g + 0(g 5 T{i 3 ) + 0(g 6 fi 4 ). 

Figure |U shows that dimensional reduction (curves labeled g 2 , g 3 , and g 4 ) ceases to be 
applicable when T < 0.2m D . At this point, non-Fermi-liquid effects take over that can 
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Figure 4. Thermal contribution to 
the interaction pressure 5p SIS db func- 
tion of T/m^ =0 for fixed chemical po- 
tential fi and coupling g = 0.5. When 
two lines of the same type run close to 
each other, they differ by changing the 
renormalization scale Ams = /i ... 4/z. 



Figure 5. The dividing line between the 
regime of dimensional reduction (DR) and that 
of non-Fermi-liquid behavior (NFL) for Nf = 3 
(full lines) and Nf = 2 (dashed lines), in com- 
parison with the weak-coupling result for the 
critical temperature of color superconductivity 
(CSC) when extrapolated to large coupling. 



only be described by HDL/HTL resummation or the new piy. In Fig. this separation 
line T « 0.2mD is plotted along with the weak-coupling result for the critical temperature 
of color superconductivity [ SJ. While at larger couplings significant non-perturbative 
modifications can be expected, at weaker couplings one observes a clear separation of 
regimes with qualitatively different weak-coupling descriptions. 
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